An efficient edge based data structure for the compressible Reynolds-averaged Navier–Stokes equations on hybrid unstructured meshes

An efficient edge based data structure has been developed in order to implement an unstructured vertex based finite volume algorithm for the Reynolds-averaged Navier–Stokes equations on hybrid meshes. In the present approach, the data structure is tailored to meet the requirements of the vertex based algorithm by considering data access patterns and cache efficiency. The required data are packed and allocated in a way that they are close to each other in the physical memory. Therefore, the proposed data structure increases cache performance and improves computation time. As a result, the explicit flow solver indicates a significant speed up compared to other open-source solvers in terms of CPU time. A fully implicit version has also been implemented based on the PETSc library in order to improve the robustness of the algorithm. The resulting algebraic equations due to the compressible Navier–Stokes and the one equation Spalart–Allmaras turbulence equations are solved in a monolithic manner using the restricted additive Schwarz preconditioner combined with the FGMRES Krylov subspace algorithm. In order to further improve the computational accuracy, the multiscale metric based anisotropic mesh refinement library PyAMG is used for mesh adaptation. The numerical algorithm is validated for the classical benchmark problems such as the transonic turbulent flow around a supercritical RAE2822 airfoil and DLR-F6 wing-body-nacelle-pylon configuration. The efficiency of the data structure is demonstrated by achieving up to an order of magnitude speed up in CPU times.     

声学超表面的非对称声分束特性研究

非对称声分束超表面是由人工微单元结构按照特定序列构建的二维平面结构,可将垂直入射的声波分成两束传播方向和分束比自由调控的透射波,在声功能器件设计及声通信领域具有广泛的应用前景。本文系统研究了一种实现非对称声分束的设计理论和实现方法,基于局域声功率守恒条件研究了声分束器的设计理论、阻抗矩阵分布、法向声强分布、声压场分布等。利用遗传算法对四串联共振腔结构进行参数优化实现了声分束器所需的阻抗矩阵分布,声压场分布表明声波入射到声分束器后在入射侧激发出两列传播方向相反且幅值和衰减系数均相同的表面波,实现了入射侧与透射侧的局域声功率相互匹配。声波经过声分束器后被分为两束透射波,两束透射波的折射角和透射系数与理论值十分吻合,证明了设计理论及实现方法的正确性和可行性。本文的研究工作可以为新型非对称声分束结构设计提供理论参考、设计方法和技术支持,并促进其在工程领域的实际应用。     

石墨基液流电池复合双极板

液流电池是一种安全性高、使用寿命长、可扩展的大规模储能系统,可以协助电网调峰储能,提高能源利用率,发展前景广阔。双极板是液流电池的重要组成部分。功能上起到了分隔、串联电池、传导电流、为电堆提供结构支撑等作用。从成本构成角度看,双极板的价格占电堆成本的比重也较大。开发高性能、低成本的双极板对加快液流电池的商业化应用具有重要意义,也是目前业界的迫切需求。虽然文献上报道了许多针对液流电池双极板开发的工作,但是目前高性能、低成本的液流电池双极板产品仍无法充分满足市场需求。本文着重介绍了石墨基复合双极板的研究现状,介绍了材料选择、工艺流程对关键性能的影响,对相关工作进行了评述,并为液流电池双极板的开发提出了建议。     

Developing Flexible Quinacridone-Derivatives-Based Photothermal Evaporaters for Solar Steam and Thermoelectric Power Generation

Solar-driven interfacial vaporization by localizing solar-thermal energy conversion to the air−water interface has attracted tremendous attention. In the process of converting solar energy into heat energy, photothermal materials play an essential role. Herein, a flexible solar-thermal material di-cyan substituted 5,12-dibutylquinacridone (DCN−4CQA)@Paper was developed by coating photothermal quinacridone derivatives on the cellulose paper. The DCN−4CQA@Paper combines desired chemical and physical properties, broadband light-absorbing, and shape-conforming abilities that render efficient photothermic vaporization. Notably, synergetic coupling of solar-steam and solar-electricity technologies by integrating DCN−4CQA@Paper and the thermoelectric devices is realized without trade-offs, highlighting the practical consideration toward more impactful solar heat exploitation. Such solar distillation and low-grade heat-to-electricity generation functions can provide potential opportunities for fresh water and electricity supply in off-grid or remote areas.     

Prediction of the drag reduction effect of pulsating pipe flow based on machine learning

Prediction of drag reduction effect caused by pulsating pipe flows is examined using machine learning. First, a large set of flow field data is obtained experimentally by measuring turbulent pipe flows with various pulsation patterns. Consequently, more than 7000 waveforms are applied, obtaining a maximum drag reduction rate and maximum energy saving rate of 38.6% and 31.4%, respectively. The results indicate that the pulsating flow effect can be characterized by the pulsation period and pressure gradient during acceleration and deceleration. Subsequently, two machine learning models are tested to predict the drag reduction rate. The results confirm that the machine learning model developed for predicting the time variation of the flow velocity and differential pressure with respect to the pump voltage can accurately predict the nonlinearity of pressure gradients. Therefore, using this model, the drag reduction effect can be estimated with high accuracy.     

Application of wake oscillators to two-dimensional vortex-induced vibrations of circular cylinders in oscillatory flows

A nonlinear time-domain simulation model for predicting two-dimensional vortex-induced vibration (VIV) of a flexibly mounted circular cylinder in planar and oscillatory flow is presented. This model is based on the utilization of van der Pol wake oscillators, being unconventional since wake oscillators have typically been applied to steady flow VIV predictions. The time-varying relative flow–cylinder velocities and accelerations are accounted for in deriving the coupled hydrodynamic lift, drag and inertia forces leading to the cylinder cross-flow and in-line oscillations. The system fluid–structure interaction equations explicitly contain the time-dependent and hybrid trigonometric terms. Depending on the Keulegan–Carpenter number (KC) incorporating the flow maximum velocity and excitation frequency, the model calibration is performed, entailing a set of empirical coefficients and expressions as a function of KC and mass ratio. Parametric investigations in cases of varying KC, reduced flow velocity, cylinder-to-flow frequency ratio and mass ratio are carried out, capturing some qualitative features of oscillatory flow VIV and exploring the effects of system parameters on response prediction characteristics. The model dependence of hydrodynamic coefficients on the Reynolds number is studied. Discrepancies and limitations versus advantages of the present model with different feasible solution scenarios are illuminated to inform the implementation of wake oscillators as a computationally efficient prediction model for VIV in oscillatory flows.     

Fluid‐structure coupling of linear elastic model with compressible flow models

Cavitation erosion is caused in solids exposed to strong pressure waves developing in an adjacent fluid field. The knowledge of the transient distribution of stresses in the solid is important to understand the cause of damaging by comparisons with breaking points of the material. The modeling of this problem requires the coupling of the models for the fluid and the solid. For this purpose, we use a strategy based on the solution of coupled Riemann problems that has been originally developed for the coupling of 2 fluids. This concept is exemplified for the coupling of a linear elastic structure with an ideal gas. The coupling procedure relies on the solution of a nonlinear equation. Existence and uniqueness of the solution is proven. The coupling conditions are validated by means of quasi‐1D problems for which an explicit solution can be determined. For a more realistic scenario, a 2D application is considered where in a compressible single fluid, a hot gas bubble at low pressure collapses in a cold gas at high pressure near an adjacent structure.     

An approach of dynamic mesh adaptation for simulating 3‐dimensional unsteady moving‐immersed‐boundary flows

In this paper, we present an approach of dynamic mesh adaptation for simulating complex 3‐dimensional incompressible moving‐boundary flows by immersed boundary methods. Tetrahedral meshes are adapted by a hierarchical refining/coarsening algorithm. Regular refinement is accomplished by dividing 1 tetrahedron into 8 subcells, and irregular refinement is only for eliminating the hanging points. Merging the 8 subcells obtained by regular refinement, the mesh is coarsened. With hierarchical refining/coarsening, mesh adaptivity can be achieved by adjusting the mesh only 1 time for each adaptation period. The level difference between 2 neighboring cells never exceeds 1, and the geometrical quality of mesh does not degrade as the level of adaptive mesh increases. A predictor‐corrector scheme is introduced to eliminate the phase lag between adapted mesh and unsteady solution. The error caused by each solution transferring from the old mesh to the new adapted one is small because most of the nodes on the 2 meshes are coincident. An immersed boundary method named local domain‐free discretization is employed to solve the flow equations. Several numerical experiments have been conducted for 3‐dimensional incompressible moving‐boundary flows. By using the present approach, the number of mesh nodes is reduced greatly while the accuracy of solution can be preserved.